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Abstract 

In  [1]  we  addressed  the  problem  of  existence  of  superconvergence  points  by  a 
computn-based  proof  and  we  gave  a  detsiled  study  of  the  superconvergence  points 
for  the  components  of  the  gradient  in  finite  element  solutions  for  Laplace’s  and 
Poisson’s  equations.  Here  we  employ  the  same  approach  to  study  the  supercon¬ 
vergence  for  the  gradient  of  the  displacement,  the  strain  and  the  stress  for  finite 
element  solutions  of  the  equations  of  plane  elasticity.  We  give  the  superconver¬ 
gence  points  for  the  components  of  the  gradient  of  the  displacement,  the  strain 
and  stress  for  meshes  of  triangles  and  squares  of  degree  p,  1  <  p  <  4.  For  the 
meshes  of  triangles  we  investigated  the  effect  of  the  topology  of  the  mesh  by  con¬ 
sidering  four  mesh-pattems  which  typically  occur  in  practical  meshes,  while  in  the 
case  of  square  elements  we  studied  the  effect  of  the  element-type  (tensor-product, 
serendipity  or  other). 
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1  Introduction 


In  [1]  we  introduced  a  computer- baaed  approach  for  finding  the  superconver¬ 
gence  points  for  the  derivatives  in  finite  element  approximations  of  Laplace’s  and 
Poisson’s  equations.  We  proved  a  mathematical  theorem  which  states  that  a  super¬ 
convergence  point  exists  if  and  only  if  it  can  be  determined  by  a  numerical  method- 
ol<^.  We  employed  the  numerical  methodology  to  find  superconvergence  points 
for  the  gradient  in  finite  element  solutions  of  Laplace’s  and  Poisson’s  equation  for 
meshes  of  triangular  or  square  elements  of  polynomial  degree  p,  for  1  <  p  <  7. 
The  conclusions  of  the  study  in  [1]  may  be  summarized  as  follows: 

1.  The  computer-based  methodology  takes  directly  into  account  the  topology 
of  the  mesh,  the  element  polynomial  spaces  and  the  type  of  the  differential 
equation  (homogeneous  or  non-homogeneous). 

2.  For  solutions  of  Laplace’s  equation  (i.e.  the  homogeneous  equation)  the  su¬ 
perconvergence  points  for  the  derivatives  always  exist  for  any  mesh-pattem 
and  type  of  elements. 

3.  For  solutions  of  Poisson’s  equation  (i.e.  the  non-homogeneous  equation)  the 
superconvergence  points  may  not  exist  depending  on  the  mesh-pattem  and 
the  element-type. 

In  this  paper  we  will  employ  an  extension  of  the  computer-based  methodology 
of  [1]  to  address  the  problem  of  existence  of  superconvergence  points  for  the  gradi¬ 
ent  of  the  displacement,  strain  and  stress  in  finite  element  approximations  of  the 
equations  of  plane  elasticity.  We  note  that  the  majority  of  the  classical  studies  on 
superconvergence  (see  [2-29]  and  the  citations  in  these  papers)  deal  primarily  with 
Poisson’s  equation  (with  a  few  exceptions;  see  for  example  [18],  [19],  [20]  and  [24] 
which  address  the  problem  of  plane  elasticity). 

The  majority  of  practical  computations  in  plane  elasticity  employ  elements  of 
degree  p,  with  1  <  p  <  4  and  most  often  p  =  2,  and  involve  the  homogeneous  case 
(i.e.  the  body-force  vanishes  identically).  Therefore  it  is  important  to  investigate  in 
detail  the  superconvergence  for  the  components  of  the  gradient  of  the  displacement, 
strain  and  stress  for  these  cases  (i.e.  for  1  <  p  <  4  and  zero  body-force).  Similarly 
as  in  [1]  we  will  study  the  superconvergence  for  the  standard  displacement  finite 
element  method  and  meshes  of  triangular  and  square  elements. 

There  are  two  types  of  superconvergence  (see  [12]),  namely: 

(i).  Direct  superconvergence:  By  this  we  mean  the  superconvergence  of  pointwise 
values  of  quantities  computed  directly  from  the  finite  element  solution. 
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(ii).  Suptrconvtrgenct  via  averaging:  The  superconvergence  of  poet-processed  val¬ 
ues  of  quantities  which  are  obtained  from  the  finite  element  solution  by  em¬ 
ploying  a  local  averaging. 

Here  (as  in  [1])  we  will  address  only  the  case  of  direct  superconvergence.  The  same 
methodologies  can  be  employed  to  study  superconvergence  via  averaging;  we  will 
present  sudi  a  study  in  a  future  paper.  Further,  to  limit  the  length  of  the  paper 
we  show  only  some  illustrative  results.  Additional  results  on  the  superconvergence 
of  various  quantities  defined  in  terms  of  the  stresses,  the  strains  etc.  can  be  easily 
obtained  using  the  approach  of  the  paper. 

Following  this  Introduction,  we  introduce  notations  for  the  model  problem  of 
plane  elasticity,  its  finite  element  approximation  and  the  types  of  meshes  employed, 
we  outline  the  theoretical  setting,  we  describe  the  numerical  methodology  for  find¬ 
ing  the  superconvergence  points  and  we  report  the  results  of  the  numerical  study 
and  the  conclusions. 

2  Preliminaries 

We  shall  consider  the  vector-valued  boundary-value  problem 


(2.1a) 

in 

(2.16) 

u,  =0 

on 

Td 

(2.1c) 

=  <. 

on 

r^v 

where  t  =  1,  2. 

Here  Q  C  is  a  bounded  domain  with  boundary  dfl  =  Pp  U  Fjy; 
n  :=  (njyn])  is  the  outward  pointing  unit-normal  on  Pj^; 

/j ,  I  =  1,  2  are  the  components  of  the  load-vector  {body-force); 

fj,  t  =  1, 2  are  the  components  of  the  normal-flux  vector  {traction)  applied  on  P 

P^  =  0,  Pp  n  Pjy  =  0;  tt  =  {u^fU^)  is  the  solution-vector  {displacement); 


(2.2a) 


,  X  1  fdui 

2\dxj'^dxi)^ 


i,j  =  1,2 
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(2.26) 


t,j  =  1,2 


kMl 

are  the  components  of  the  strain  and  the  stress,  respectively; 

°>iu>  *>  it  6,  /  =  1,  2,  are  the  material-coe£Bcients  {elastic  constants)  which  satisfy 


(2.3a) 


<*y  W  =  ^jilk  =  *1  i,  /  =  1 , 2 


3  3 

(2.36)  ^ijkl^ij  ^kl  —  ^  ®  ^  ® ^  ^ij  ~  ^ji 

(Conditions  (2.3a),  (2.3b)  are  satisfied  for  linear  anisotropic  elasticity;  in  the  case 
of  isotropic  plane  elasticity  a^^y  =  p(6<j6w+6a6*,)  where  6^^  is  Kronecker’s 

delta  and  A,  p  are  Lame’s  constants.) 

Let  us  now  cast  the  model  problem  in  variational  form.  We  will  use  the  nota¬ 
tions 


(2.4a) 


H*  :=  I  e  =  (ui,U2)  *  ^  } 


(2.46) 


==  {  ®  =  (''i»*'3)  1  V.  €  if*(n),  V,  =0  on  } 


(2.4c)  \\v\\ija  :=  ||w.  ||?,n)  ,  |v|,,n  :=  l?.n) 

with  ||v,-  IIj  n  (resp.  |v,-  |i^)  being  the  usual  H^{0)  Sobolev  norm  (resp.  seminorm). 
The  variational  form  of  the  boundary- value  problem  (2.1)  is  now  posed  as: 


Find  u  €  ^ch  that 

*  D 


2  2 

(2.5a)  Bn(t*, «)  =  /  S  /.  «<  +  /  E  ^  €  Hj. 

tal  •'»  ir  tal 


where  the  bilinear  form  x 


R  is  defined  by 


(2.56) 


The  energy-norm  over  any  subdomain  5  C  (2  is  defined  by 


(2.5c)  |||v|||5  :=  /Bs{v,v) 

where  Bg{UjV)  has  the  obvious  meaning. 

Let  T  =  {T^}  be  a  family  of  meshes  of  triangles  or  quadrilaterals  with  straight 
edges.  It  is  assumed  that  the  family  is  regular^  namely:  For  the  triangles  the 
minimal  angle  of  all  the  triangles  is  bounded  below  by  a  positive  constant,  the 
same  for  all  the  meshes.  For  the  meshes  of  quadrilaterals  it  is  assumed  that  the 
mesh  can  be  mapped  to  a  mesh  of  squares  by  a  sufficiently  smooth  transformation 
and  hence  it  is  sufficient  to  study  the  superconvergence  for  a  mesh  of  squares  (see 
[9],  [10]  for  details).  Let  us  introduce  the  finite-element  spaces 

(2.6)  s;:={«€H*|o,.oy,^€5’{f),  i  =  l,2,  i  =  l . 

where  is  the  miq>ping  function  for  the  ibth  finiteelement  which  maps  either  a 
standard  triangular  element  (using  a  linear  transformation)  or  a  standard  quadri¬ 
lateral  element  (using  a  bilinear  transformation)  onto  the  ibth  finite  element,  f 
denotes  a  standard  element,  Ai{T|^)  is  the  number  of  elements  in  the  mesh  T|^y 
denotes  the  element-space  over  f . 

As  in  [1]  we  will  consider  the  following  choices  for  the  element-space  ^(r): 
a.  Complete  polynomial  space  up  to  degru  p. 

For  the  triangular  elements  we  let  5p(f)  =  Vp{f)  where 

(2.7)  ^,(#):={>|p(i„i,)= 

'  * 


For  the  square  elements  we  consider  the  following  choices  for  the  definition  of 
the  polynomial  space  5^(^)  (see  also  [32]). 
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h.  Tensor-produet  (bi-p)  polynomial  spaix  of  degree  p. 


(2.8)  i<">(#):={p|P(i„i,)= 


c.  Serendipity  (trunc)  polynomial  space  of  degree  p. 

(2.9)  >(f )  :=  (P|  >(i.,  i,)  =  E  “«  ij  +  *>  +  “..p  * .  ^  } 

V  *4  •' 

d.  Intermediate  polynomial  space  of  degree  p. 


p-i 


(2.10)  5"(f);=|>|>(f„ij)=  + 


M 

®s*+>s» 


We  let 


(2.11) 

The  finite  element  solution  of  the  elasticity  problem  satisfies: 
Find  U|^  € 


(2.12) 


3  3 

5n(«fc.  =  /,  E  fi  "w  +Ltl9iVHi  V  6 


The  error  is  e^\=u  —  u,^. 

3  Definition  of  the  superconvergence  quantities 

Let  €  T  be  a  finite  element  grid,  the  corresponding  finite 

elemoit  solution  and  let  t  €  be  any  element.  Let  F(u)  be  the  solution  quantity 

of  interest,  for  example  F{u)  =  ^  or  €<,(«)  or  i,  j  =  1,  2.  We  are 

interested  in  the  values  of  the  relative  error  in  F’(u)  at  points  €  r, 
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if  «(tt  -  ttfc)  9^  0 

if  'i’(u  —  =  0 


where 


(3.U) 

Let  us  assume  that  r  is  an  element  of  fixed  geometry  (but  not  of  fixed  size) 
which  appears  in  all  the  meshes  of  the  sequence  T.  If  there  exists  a  point 
which  is  fixed  with  respect  to  the  element  r,  such  that 

(3.2)  (im  e(*^;  F;  t»,  «*,/»,  t)  <  ,  0  <  »;  <  100 

we  will  say  that  is  a  u  —  ri%~auperconvtrgtnct  point  in  element  t.  Further,  we 
will  call  x^nU— i;%-8uperconvergcnce  point  if  it  is  u  —  17%-superconvergence  point 
for  every  u  £  U.  Note  that  0(«;  F;  u,  and  thus  all  points  in  every 

element  r  are  100%-8uperconvergence  points.  If  there  exists  a  point  x^  such  that 

(3.2)  holds  for  a  given  solution  u  (resp.  class  of  solutions  U)  with  >7  =  0  then  x^ 
is  a  u-superconvergence  (resp.  2Y-8uperconvergence)  point  in  the  classical  sense. 

Let  us  now  define  the  following  geometrical  quantities  which  will  be  employed 
in  the  study  of  superconvergence  below.  For  a  ^ven  >7,  0  <  >7  <  100  we  define 

1.  rj%-contour  of  F{u)  in  the  element  r  of  the  mesh  for  the  exact  solution  u: 

(3.3)  :=  {*  €  r  I  ©(*;  F;  u,  u*,  h,  r)  =  | 

2.  Superconvergence  points  of  F(u)  in  the  element  r  of  the  mesh  Ti^  for  the  class 
of  exact  solutions  U: 

(3-4)  n  n)  :=  n  <%)(»:’■.  3a) 

«€</ 
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S.  ri%-band  of  F{u)  in  the  dement  r  of  the  mesh  for  the  exact  solution  u: 


(3.5)  T,  Th)  :=  I*  €  t  I  e(*;  F;  u,  *’’’)<  } 

4.  tf%-superconvergence  regions  of  F{u)  in  the  element  r  of  the  mesh  for  the 
doss  of  exact  solutions  U\ 

(3.6)  r,  n)  :=  fl  «??.)(«:  n  Tk) 

uCU 

Remark  S.l.  We  could  be  interested  in  the  points  which  are  superconvergent  simul¬ 
taneously  for  several  functionals  e.g.  all  the  stress  components.  We  can  formalize 
this  by  assuming  a  vector  functional;  the  meaning  of  the  t7%-superconvergence 
point  for  the  vector  functional  is  obvious. 

4  The  class  of  locally  periodic  meshes 

In  this  paper  we  will  study  the  superconvergence  of  finite  element  solutions  for 
the  equations  of  plane  elasticity  for  a  special  class  of  locally  periodic  meshes  which 
are  d^ned  as  follows.  Let  0  <  H  <  H°  =  (x®,  x®)  6  , 

(4.1)  5(*°,  H)  :=  I*  =  (xi,  X2)|  |xf  -  x^\  <  H ,  i  =  1, 2| 

and  assume  /f®  is  sufficiently  small  such  that  5(x®,  /f®)  C  H.  Further,  let  7  be  a 
set  of  multi-indices  {i,j) ,  =  (xj'’'^  Xj  ■'^)  €  11  and 

(4.2)  c(*(‘-'J,/t):=5(*(‘^,/i)c5(*®,/f),  (i,i)  €  7 

be  the  set  of  the  h-cells  (or  cells)  which  cover  exactly  S{x°,H)  i.e. 


(4.3a)  (J  c  (*(••'),  h)  =  5(*®,^) 

(•j)€7 


(4.36)  c(*^“•'‘^6)f)c(*<‘’•"^A)  =  0  for  (t'l.ji)  #  (tj.ia) 
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We  will  refer  to  S(x^f  H)  as  the  subdomain  of  periodicity  of  the  mesh  centered  at 
We  will  denote  by 

(4.4)  c  :=  5(0,1)  :=  |(ic,,i3)|  jx,|  <  1 ,  |x2|  <  l} 

the  unit-  {master-)  cell  c,  the  h-cell  is  an  h-scaled  and  translated  master-cell. 

Let  T  be  a  mesh  of  triangles  or  squares  on  the  master-cell  (the  master-mesh) 
and  be  the  mesh  on  c{x^''^\h)  which  is  the  scaled  and  translated  image 
of  T.  We  will  consider  the  family  T  of  locally  periodic  meshes.  Let  Tf^  £  T  and 
r/^(*®,  H)  be  the  restriction  of  on  5(*, H)  and  the  restriction  of  T^^{x°,  H) 
on  c(*^*’^^fc).  We  assume  that  (i,j)  €  7  i.e.  Ti^{x°,H)  is  made  by 

the  periodic  repetition  of  the  h-scaled  master  mesh. 

The  type  of  meshes  imder  consideration  is  depicted  in  Fig.  1,  where  the  periodic- 
mesh  subdomaun  S{x'‘^,H)  is  shown  with  thick  perigram.  Outside  the  subdomain 
5(«°,  H)  the  mesh  is  arbitrary;  it  could  have  curved  elements,  refinements,  etc. 


5  Outline  of  the  theoretical  setting 

Let  Q  be  a  vector- valued  function  with  components  which  are  polynomials  of 
degree  (p  +  1)  defined  over  the  master-cell  c  and  let  T  be  the  master-mesh.  Then 
denote 


(5.1) 

where  is  the  interpolant  of  degree  p  of  the  function  Q  defined  over  the  muster, 

mesh  T  (for  which  h  ^  1).  Any  vector- valued  function  with  components  which 
are  polynomiab  of  degree  p  on  an  element  belongs  to  5i(rj^)  and  hence  any 
polynomial  of  degree  p  on  S{x°,H)  belongs  to  Sf(S(x°,If)).  It  follows  that  p 
defined  in  (5.1)  is  c-periodic;  this  can  be  shown  exactly  as  in  [40].  We  have 


(5.2o) 

p(l,X3)=  p(-l,X3),  |X2l<l 

(5.26) 

p(xi,l)=  p(i„-l),  |x,|<l 

Let 
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(5.3) 


Hpbr(c)  :=  |tt  €  H'(c)  tt  satisfies  (5.2)  | 


and 


(5.4)  S;j.er(S)  :=  {«  €  fftEl>(S)|  »|,  6  St(f) ,  i  =  1,2,  V  f  €  f } 
F\irther  let  €  ^pbr(c)  such  that 

(5.5a)  B^(x^«)=  Bg(p,v)  V  r  6  SJ,per(c) 

and 

(5.5i)  jf  (p  -  i')  =  0 

Note  that  the  function  exists  and  is  uniquely  detennined  (we  will  compute  it 
numerically  in  the  examples).  Let  us  also  define  fjf  €  H*(c)  by 

(5.6)  ff> p  —  z*^  Q  —  w  where  w  :s=  +  x'’ 

Let  V’a  €  h))  be  the  function  defined  above,  scaled  and  translated 

onto  the  cell  c(*bJ)^ftj  of  the  mesh  in  S{x°,H)  i.e. 

(5.7)  :=  ,  ^(*)=  i  =  l,2, 

where  *  =  t(*  ~  *  €  c(ajb*>),h).  It  is  easy  to  see  that  V’a  con  be 

h 

periodically  extended  over  5(sb°,  ). 

In  [1]  we  proved  the  following  theorem  for  Poisson’s  equation  based  on  the 
theory  of  interior  estimates  (see  [33]-[39]): 

Theorem  1.  (Poisson’s  equation;  see  [1])  Let  <  H  <  and  assume  that  the 
following  assumptions  hold  with 


/9=  p+1  -c, 


1 

6(6p+l) 


(5.8) 


6p+  1 
6p  ’ 
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c  = 


Assume  th&t  the  exact  solution  u  satisfies 


(5.9a) 


<  K  <  oOf  0<|a|<p  +  2 


where  a  :=  (a,, a,),  D^u  :=  ,  I<»1  :=  <*1  +  «a. 


(5.96)  ^  a*  >  0  where  :=  (Z)®u)(x®) 

Further  assume  that  the  mesh  T|^  is  such  that 

(5.10)  lkfcllt>(5(xO.H,))  ^  Ch^Hi ,  with  ^  >  (p  +  1)  —  c 
Moreover  assume  that  the  meshes  in  5(x®,  H)  are  such  that 

(5.11)  CiHt<h<C^Ht 
Then  for  any  x  €  5(*®,  Hi ) 

(5.12)  (,)|  =  1^  (,)|  +  aC5'«« 
with  |A|  <  1  and  C  independent  of  h. 

Remark  5.1.  The  theorem  assumes  that  the  mesh  is  periodic  in  a  small  subdo¬ 
main  (i.e.  5(c®,  H))  in  the  interior  of  the  domain  and  that  the  solution  is  smooth 
in  the  neighborhood  of  the  subdomain.  Outside  the  subdomain  we  assume  nei¬ 
ther  periodicity  of  the  mesh  nor  smoothness  of  the  solution.  The  solution  may 
have  algebraic-type  singularities  at  a  finite  ntunber  of  comer  points  or  points  of 
abrupt  change  in  the  type  of  boundary-condition.  Here  it  is  only  assumed  that  the 
pollution-error  in  a  shrinking  mesh-patch  (i.e.  T|^{x'^^Hl))  in  the  interior  of  the 
subdomain  is  controlled;  this  implies  that  the  mesh  has  been  adequately  refined  in 
the  neighborhood  of  all  singular  points. 

Remark  5.2.  If  we  further  assume  that 


t«(T) 


>Ch^  *  =  1,2,  c>0 


theorem  1  implies  that:  A  point  in  the  element  r  is  a  superconvergence  point 
for  Fiv)  —  -T —  in  the  element  t  if  and  only  if  -z — (as^)  =  0.  Assumption  (5.13) 

aXi  oXi 

can  be  realized  by  imposing  additional  restrictions  on  the  values  of  the  (p  + 1)- 
derivatives  of  the  solution  at  Xq.  This  assumption  is  reasonable  because  we  are 
interested  in  a  sufficiently  large  class  of  solutions  U. 

Remark  5.S.  Under  assumption  (5.13)  Theorem  1  also  states  that:  A  point  Xf 

in  the  element  f  is  an  asymptotically  Ti%-superconvergence  point  for  F{u)  =  — , 

uX^ 

i  =  1,2,  if  and  only  if  Q(Xf;F;QyW,l,f)  <  where 


(5.14)  e(a!,;  F]  Q,  w,  1,  f)  := 


l^(V>)(»f)l 

W) 

0 


,  if  '1(0)  :=  ||0||i«(^)  0 

1  if  IIV'|lL~(f)  =  0 


Remark  S.f.  The  proof  of  theorem  1  in  [40]  was  based  on  various  interior  estimates 
for  the  error  in  finite  element  approximations  of  Poisson’s  equation,  especially  the 
results  given  in  [38]  and  [39].  It  is  very  plausible  that  analogs  of  these  results 
hold  for  finite  element  approximations  of  the  elasticity  equations  and  more  general 
elliptic-systems  because  the  main  ideas  of  the  proofs  of  these  results  carry  to  the 
general  case.  To  our  knowledge  the  precise  details  are  not  available  for  the  elasticity 
equations.  Nevertheless  we  will  assume  the  validity  of  the  analog  of  Theorem  1  for 
the  equations  of  elasticity. 


6  The  methodology  for  determining  the  super¬ 
convergence  points 

In  order  to  study  the  superconvergence  of  finite  element  solutions  in  uniform 
mesh-patches  in  the  interior  of  a  subdomain  5(x°,  /f),  we  let 

(6.1)  :=  {u  €  H'(ll)|  <  If ,  i  =  1,2,  0  <  |<»|  <  p  +  2] 

the  class  of  solutions  which  are  locally  smooth  in  5(x°,  H),  where  5(x°,  H)  denotes 
an  interior  subdomain  of  interest  in  which  the  mesh  is  locally  periodic  as  described 
above  (the  subdomain  must  be  a  finite  distance  away  from  the  boundary  and  points 
of  roughness  of  the  body-force;  see  Fig.  1).  In  the  majority  of  the  applications  one 
is  only  interested  in  the  subclass  of  solutions  in  14^  which  are  “harmonic”,  namely. 
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(6.2) 


«■*■;=  I  4(«)=0,  i  =  l,2,  in  o} 

We  may  also  assume  that  the  functions  are  “hannonic”  in  a  subdomain  which  is 
slightly  bigger  than  5(s;°,  H)  and  which  includes  S{at^y  H)  in  its  interior. 

For  a  given  locally  periodic  grid  with  corresponding  periodic  master-mesh  T, 
given  material  orthotropy  and  given  class  of  smooth  solutions  U  we  let 

(6-3)  C;=  f«|  «(i„i»)  =  L  “»«»(*>. *2)) 

denote  the  class  of  (p  +  l)-drgree  monomials  which  occur  in  all  (p  -f-  l)-degree 
Taylor-series  expansions  of  functions  frmn  U.  Here  Qi,,  k  —  1, . . . ,nd  denotes  a 
set  of  linearly  independent  monomials  which  form  a  basis  for  Q.  For  example,  let 
us  assume  that  If  is  the  class  of  smooth  solutions  given  in  (6.1);  in  this  case 
Q  is  the  2(p  -t-  2)  dimensional  space  of  vector-valued  functions  with  components 
which  are  monomials  of  dq;ree  (p+  1).  The  set  Q  which  correspondi*  the  class 
of  “harmonic”  solutions  14“^  is  the  four-dimensional  linear  space  of  “harmonic” 
monomials  of  degree  (p+ 1)  denoted  by  The  “harmonic”  basis  monomials  of 
degree  (p  +  1)  for  1  <  p  <  4,  which  were  employed  in  the  computations,  are  given 
in  the  Appendix. 

The  asymptotic  values  of  the  error  for  any  smooth  solution  u  in  the  interior  of  a 
periodic  mesh-subdomain  can  be  obtained  by  solving  the  periodic  boundary- value 
problem  (5.5),  using  the  master-mesh  T  over  the  master-cell  c,  with  data  obtained 
from  the  local  (p-h  l)-degree  Taylor-series  expansion  Q  of  the  exact  solution.  The 
asymptotic  i7%-contour8  for  a  given  solution  u  can  be  obtained  by  contouring  the 
function  F(^),  with  ^  defined  as  in  (5.6)  corresponding  to  the  local  Taylor-series 
expansion  Q  of  the  solution  «.  The  superconvergence  points  x  for  a  given  class  of 
solutions  If  satisfy 


(6.4)  F{^i){x)  =0,  1  <  i  <  nd 

Therefore  «  is  a  superconvergence  point  if  and  only  if  the  zero-contours  of  F(^,  ) 
intersect  at  «  for  1  <  >  <  nd.  Here  :=  p,  —  which  is  obtained  from  (5.6) 
for  Pi  s=  Qi  —  {Qi  )f^  where  Q,-  is  the  i-th  basis  monomial  of  the  nd-dimensional 
monomial  space  Q  corresponding  to  the  class  U.  We  also  let  to,  :=  (Q,-  )f^  -f  . 

The  asymptotic  17%-superconvergence  regions  for  a  class  of  solutions  U  can  be 
determined  by  using  numerical  optimization.  In  particular,  let  us  consider  the 
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unifOTm  subdivision  of  the  element  r  into  subtriangles  with  vertices  at  the  set  of 
pmnts  E  :=  define  the  function 


(6.5)  F;  Q,  {&,  1,  f)  :=  max  - -  100 , 

‘  I  max  |Ea.  W.)(e,)l 

at  the  points  in  E.  The  function  will  be  defined  for 

any  point  «  €  f  by  using  linear  interpolation  in  the  subtriangles.  The  asymptotic 
17%-superconvergence  regions  in  the  element  r  can  be  approximated  using  the  level- 
sets  of  the  functions  F;  1,  f)  i.e. 

(6.6)  f ,  f )  »  {i  e  f  I  F;  C,  {&,  )S„  1 ,  ? )  <  .)%} 

We  will  call  the  above  approach  the  direct  approach.  It  is  also  possible  to  use  a 
simplified  approach  which  avoids  the  use  of  numerical  optimization  at  every  point. 
Let  us  define  (see  also  [30,  31]) 


(6.7) 

where 


Zs  :=  min 


-  .si _ 


(gof)i 

•si  . 


The  quantity  Zs  can  be  computed  using  numerical  optimization.  Then 


(6.9)  *J|.,(C;  r,  f )  :=  {i  6  r  |  ^.,(*;  F;  C.  {*i  1,  f )  <  } 

are  the  sqiproximate  regions  of  q%-superconvergence  for  the  class  of  solutions  Q. 
Remark  6.1.  Note  that 


(6.10) 


^n.)(C;^7r)c^^^^,((2;f,r). 
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Therefore  the  simplified  approach  results  to  a  conservative  estimate  for  the  17%- 
superconvergence  regions. 

Remark  6.S.  The  functions  defined  in  (6.5),  (6.7)  depend  on  the  set  of  points  E. 
To  ensure  good  accuracy  in  the  approximation  of  the  t7%>superconvergence  regions 
a  sufficient  number  of  points  must  be  employed. 

7  Numerical  study  of  superconvergence  for  pe¬ 
riodic  meshes  of  triangles  and  squares 

We  will  now  use  the  methodology  of  the  previous  Section  to  find  the  super¬ 
convergence  points  or  the  17%-superconvergence  regions  for  the  components  of  the 
gradient  of  the  displacement,  strain  and  stress  for  finite  element  solutions  of  the 
equations  of  elasticity  in  the  interior  of  periodic  meshes  of  triangular  and  square 
elements.  In  the  numerical  examples  we  addressed  the  following  questions: 

1.  For  periodic  meshes  of  triangles  with  various  mesh-topologies,  and  elements 
of  degree  p,  where  are  the  superconvergence  points  for  the  various  solution 
quantities  for  the  class  of  "harmonic”  solutions?  Are  these  points  supercon¬ 
vergence  points  for  the  class  of  general  solutions? 

2.  For  meshes  of  squares  of  degree  p,  where  are  the  superconvergence  points  for 

the  various  quantities  for  the  tensor-product  space  ),  the  serendipity 

space  ^(r)  and  the  intermediate  space  S^{t)7 

3.  In  the  cases  that  there  are  no  superconvergence  points  (i.e.  0%-superconvergence 
points)  where  are  t7%-superconvergence  points  for  small  values  of  17%? 

We  will  answer  these  questions  using  the  computer-based  approach  of  Section  6. 

7.1  Determination  of  the  superconvergence  points  for  the 
periodic  meshes  of  triangles 

The  majority  of  the  results  for  the  superconvergence  points  for  the  triangular 
elements  in  the  literature  are  given  exclusively  for  the  Regular  pattern  (which 
is  also  known  as  the  three-directional  mesh  and  is  shown  in  Fig.  2a),  for  linear 
and  quadratic  elements  and  for  the  Poisson's  equation.  In  [1]  we  determined  the 
superconvergence  points  for  Laplace’s  and  Poisson’s  equations,  for  all  the  mesh 
patterns  shown  in  Fig.  2  and  elements  of  degree  p,  1  <  p  <  7.  Here,  we  employed 
the  numerical  methodolo^of  Section  6  to  find  the  superconvergence  points,  for  the 

same  mesh-pattems,  for  for  the  class  of  "harmonic”  solutions  of  the  equations 

of  plane-elasticity  and  1  <  p  <  4. 
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In  Fig.  3  (reap.  Fig.  4)  we  give  examples  of  how  the  superconvergence  points 
for  the  class  of  ‘‘harmonic”  solutions,  for  p  =  1, 2,  3,  were  obtained  for  the  R»ular_ 
(resp.  Criss-Cross)  pattern  from  the  intersection  of  the  zero-contours  (Qji  ^)t 

of  the  error  functions  0,-,  t  =  which  correspond  to  the  basis  "harmonic” 

monomials  of  degree  (p  -h  1).  FVom  the  numerical  results  we  observe  that: 

(i)  In  the  Regular,  Chevron  and  Criss-Cross  patterns  for  linear  or  cubic  elements 

(p—  1  or  3),  there  exists  one  superconvergence  point  for  in  the  elements 

with  an  edge  parallel  to  the  x^-axis.  This  point  is  located  at  the  midside  of 
the  edge  parallel  to  the  X|-axis. 

(ii)  In  the  Regular,  Criss-Cross  patterns  for  quadratic  elements  (p  =  2)  there  are 
two  superconvergence  points  for  in  the  elements  with  an  edge  parallel  to 

(/Xj 

the  X|-axis.  These  points  are  located  at  the  Gauss-points  of  the  edge  parallel 
to  the  X|-axis. 

(iii)  In  the  Chevron  pattern,  the  Union-Jack  pattern  and  element  Tj  of  the  Criss- 
Cross  pattern  (shown  in  Fig.  2d)  for  quadratic  elements  (p  =  2)  there  are  no 

dui 

superconvergence  points  for 

i/Xj 

(iv)  The  superconvergence  points  for  the  components  of  the  gradient  of  the  dis¬ 
placement  for  the  class  of  general  solutions  coincide  with  the  superconver¬ 
gence  points  for  the  class  of  "harmonic”  solutions. 

(v)  There  are  no  superconvergence  points  for  the  stress  in  any  of  the  mesh- 
pattems  except  for  the  special  case  that  the  Poisson’s  ratio  is  equal  to  zero. 

In  this  case  the  normal  stress  components  an,  a^i  are  superconvergent  at 

the  superconvergence  points  for  respectively  (if  such  points  exist 

in  the  mesh-pattem). 

(vi)  There  are  no  superconvergence  points  for  the  shear-stress  in  any  of  the  mesh- 
pattems. 

(vii)  For  p  —  4  there  are  no  superconvergence  points  in  any  of  the  mesh-pattems 
for  any  of  the  quantities  and  for  all  values  of  Poisson’s  ratio. 
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7.2  f7%>superconvergeiice  regions  for  the  components  of 
the  gradient  and  stress  for  the  periodic  meshes  of 
triangles 


For  p  >  3  there  are  very  few  (if  any)  superconvergence  points  for  any  of 
the  solution-quantities  for  the  problem  of  plane  isotropic  elasticity.  In  these  cases 
suitable  sampling  points  (i.e.  points  where  the  error  in  the  solution-quantity  is 
small,  asymptotically,  with  respect  to  the  error  in  other  points  in  the  element) 
can  be  determined  from  the  f7%-superconvergence  re^ons.  In  the  Figs.  5-7  below 
we  give  the  regions  TCSt  (2“^;^;7),  for  the  Regular,  Chevron  and  Criss-Cross 

patterns  for  p  =  1,  2,  3  and  Poisson’s  ratio  —  0.3.  (It  should  be  noted  that  the 
17%-scale  employed  varies  in  the  Figures.)  FVom  the  numerical  results  we  observe 
that: 


(i)  The  regions  ^^(C^;t;T)  exist  for  small  q  for  almost  all  the  cases  (note 

however  that  rfk  =  100%  everywhere  in  the  element  Tj  of  the  Criss-Cross 
pattern  for  linear  elements,  as  shown  in  Fig.  7a). 


(ii)  For  p  =  2,  in  the  Chevron  pattern  there  exist  q%-superconvergence  regions 
with  minimum  q  fts  34%  (as  shown  in  Fig.  6b). 

(iii)  For  p  =  2,  3,  in  the  element  r,  of  the  Criss-Cross  pattern  there  exist  q%- 
superconvergence  regions  of  significant  size  (although  q  may  be  relatively 
large),  as  shown  in  Figs.  7b,  7c,  respectively. 

FVom  the  engineering  point  of  view  one  is  mostly  interested  in  determining 
optimal  sampling  points  for  the  stress-components.  Except  for  the  special  case 
of  zero  Poisson’s  ratio  there  are  no  superconvergence  points  for  any  of  the  stress- 
components  in  any  of  the  mesh-pattems  and  for  elements  of  any  degree  p.  Here 
we  show  that  for  elements  of  degree  p  >  2  there  exist  q%-superconvergence  points 
for  the  stress-components  in  all  the  mesh-pattems  for  relatively  small  values  of 
q;  these  points  may  be  employed  as  sampling-points  for  the  corresponding  stress- 
components.  In  Figs.  8, 9, 10  we  give  the  q%-superconvergence  regions  for  Ou  and 
<Ti2,  for  p  =  2  and  3,  for  the  Regular,  Chevron  zmd  Criss-Cross  patterns.  (The 
regions  for  the  Union- Jack  are  similar  and  will  not  be  given  here;  note  that  the 
Union-Jack  pattern  is  obtained  by  a  45*-rotation  of  the  Criss-Cross  pattern.)  We 
did  not  give  the  regions  for  p  =  1  because  the  minima]  values  of  q%  for  which 
the  regions  exist  are  close  to  100%  for  all  the  mesh-pattems  (in  other  words  for 
meshes  of  linear  triangles  there  is  not  a  preferable  set  of  sampling  points  for  the 
stress-components,  unless  the  Poisson’s  ratio  is  equal  to  zero). 

The  minimal  value  of  q%  in  the  q%-superconvergence  regions  shown  in  Figs.  8- 
10  are  given  in  Table  1.  (We  did  not  report  the  points  where  these  minima  occur 
but  these  can  be  easily  found.)  In  summary  we  observed  that: 
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(i)  The  regions  T)  for  F(u)  —  Oj,  or  o,,  exist  (resp.  do  not  exist) 

for  p  =  2,  3  (re^.  for  p  =  1)  for  relatively  small  values  of  17%. 

(ii)  The  minimal  values  of  17%  and  their  locations  depend  on  the  mesh-pattem, 
the  degree  p  of  the  elements,  the  stress-component  and  the  orientation  of  the 
coordinate-axes  with  respect  to  the  mesh-pattem. 


7.3  Superconvergence  points  for  periodic  meshes  of  squares 

We  also  used  the  computer-  based  approach  to  determine  the  superconver¬ 
gence  points  for  meshes  of  square  elements.  Here  we  investigated  the  effect  of  the 
choice  of  the  finite-element  space  (tensor-product,  serendipity  and  intermediate 
element-space)  on  the  superconvergence  points  for  the  components  of  the  gradient 
of  displacement,  strain  and  stress.  In  the  cases  where  there  are  no  superconver¬ 
gence  points  (for  example,  in  the  quartic  serendipity  element)  we  reported  the 
element-coordinates  of  f7%-superconvergence  points  for  minimal  values  of  17%. 

o.  Tensor-product  and  initrmtdiait  family 

dn  dti 

The  superconvergence  points  for  and  are  located  on  the  Gauss-lines 

which  are  parallel  to  the  x^-axis  and  are  intersecting  the  Xj-axis  at  the  Gauss- 
Legendre  points  of  degree  p.  The  (p  x  p)  Gauss-Legendre  points  are  superconver¬ 
gence  points  for  all  the  components  of  strain  and  stress  for  all  Poisson’s  ratios  for 
both  the  “harmonic”  and  general  class  of  solutions. 

b.  Senndipity  family 

For  p  =  1  and  p  =  2,  the  superconvergence  points  for  all  the  quantities  are 
exactly  the  same  as  the  corresponding  sunerconvergence  points  for  the  tensor- 
product  and  the  intermediate  family  i.e.  the  p  x  p  Gauss-Legendre  points  in  the 
element. 

In  the  cubic  serendipity  element  (p  =  3)  there  exist  four  superconvergence 
points  and  one  superconvergence  line  for  the  components  of  the  gradient  of  the 
solution  for  the  class  of  “harmonic”  solutions.  The  superconvergence  points  for 
dU| 

— —  are  given  in  Table  2  and  are  shown  in  Fig.  11a.  The  four  superconvergence 
ax, 

dtA 

points  and  the  superconvergence  line  for  are  shown  in  Fig.  11b.  For  Cjs  (and 

(/X2 

(T12)  and  Oji  (and  there  is  only  one  superconvergence  point,  at  the  center  of 

the  element,  as  shown  in  Figs.  11c  and  lid,  respectively.  All  points  given  above  are 

also  superconvergence  points  for  the  class  of  general  solutions  and  all  admissible 

values  of  Poisson’s  ratio.  In  the  special  case  of  zero  value  for  the  Poisson’s  ratio 

the  superconvergence  points  for  (resp.  coincide  with  the  superconvergence 

du^  .  du^ 
points  for  (resp.  -r-^). 

vX|  0X2 


For  the  cubic  serendipity  element  we  also  determined  the  17%-superconvergence 

regions  for  F{u)  =  <t„  and  (r„;  these  are  given  in 

Figs.  12a,  12b,  12c,  12d,  respectively,  for  Poisson%  ratio  equal  to  0.3.  Prom  Fig.  12c 
we  observe  that  some  of  the  points  of  the  3x3  Gauss- Legendre  product  rule  (these 
points  are  often  employed  to  sample  the  stresses  in  the  cubic  serendipity  element) 
correspond  to  17%-superconvergence  points  for  the  normal  stress  components  with 
17%  >  75%.  From  the  same  Figure  it  is  clear  that  it  is  possible  to  find  sets  of 
sampling  points  for  the  normal  stresses  with  17%  <  50%. 

In  the  quartic  serendipity  element  (p  =  4)  there  arc  no  0%-superconvergence 
points  for  any  of  the  solution  quantities.  We  found  however  that  there  exist  17%- 
superconvergence  points  and  regions  for  the  components  of  the  gradient,  strain 
and  the  stress  (for  v  =  0.3)  for  small  values  of  17,  namely  17%  <  2.5%.  In 

Figs.  13a  and  13b  we  give  the  J7%-superconvergencc  regions  (Q“^;T,!f)  and 

respectively,  for  Poisson's  ratio  equal  to  0.3.  In  Tables  3b  and 

3c  we  give  the  master-element  coordinates  of  sampling-points  for  £13  and  with 
17%  <  2.5%.  The  q%-superconvergence  regions  for  ^  25  are  shown  in 

Figs.  13c,  13d,  respectively.  We  also  determined  the  common  i7%-8uperconvergence 
regions  for  a  class  of  Poisson's  ratios  (0  <  v  <  0.35).  In  Fi^s.  14a,  14b  and  14c 

we  show  the  regions  f)  '^fM«)(C*^;^;7’)  for  F{u)  =  <y,i  and  fjj,  re- 

OKyKOJK  . 

spectively,  for  p  =  4  and  the  class  of  “harmonic”  solutions  (the  Poisson's  ratio  was 
varied  from  0  to  0.35  in  steps  of  0.05).  FVom  Fig.  14b  it  can  be  seen  that  there 
is  a  very  small  common  25%-supcrconvergence  region  for  <Tjj  for  all  the  Poisson's 
ratios  (0  <  1/  <  0.35). 


7.4  Rate  of  convergence  at  the  superconvergence  points 

0XA 

We  checked  the  rate  of  convergence  of  at  the  superconvergence  points 

(given  in  Sections  7.1  and  7.3)  in  model  computations  using  relatively  coarse 
meshes.  We  considered  the  Dirichlet  problem  with  data  consistent  with  the  exact 
solution  =  U3(x,,i3)  =  sin(TXi)sin(xz3)  (note  that  tt  ^  W“^)  in  the 

domain  ft  =  (0, 1)’  which  was  meshed  by  a  uniform  grid  of  elements  (of  triangles 
in  the  Regular-pattern  or  squares).  We  computed  the  quantity: 


where  denotes  the  superconvergence  points  in  the  elements  in  the  subdomain 
fto  =  (0.25,0.75)’.  We  computed  the  values  of  E  in  meshes  with  mesh-sizes  h  = 
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7*  A*  values  of  £  are  superconvergent  with  rate  (p  +  (Tq) 

4  8  16 

if  there  exists  (Tq  >  0  s.t.  =  constant. 

In  Table  4a  (resp.  Table  4b)  we  give  the  values  of  £  and  h~^£  (resp.  h~*E) 
computed  using  a  uniform  mesh  of  triangles  of  Regular  pattern  with  p  =  2  (resp. 
p  =  3).  We  note  that  as  the  element  size  h  is  decreased,  the  values  of 
converge  to  constants;  therefore  the  quantity  £  converges  at  the  rate  of  (p+ 1)  i.e. 
Oq  =  1.  In  Table  4c  we  give  the  values  of  £  and  h~*E  computed  using  a  uniform 
mesh  of  cubic  serendipity  squares  of  size  h.  We  computed  the  values  of  £  using  the 
points  given  in  Table  2  in  the  subdomain  Oq.  It  can  be  observed  that  the  values  of 
the  quantity  £  converge  with  rate  equal  to  4.  Hence  the  points  given  in  Sections 
7.1  and  7.3  are  superconvergence  points  with  rate  equal  to  (p  +  1). 

We  also  checked  the  value  of  relative  error  at  the  t7%-superconvergence  points 
for  the  quartic  serendipity  element  which  are  given  in  Table  3a.  We  considered 
the  Dirichlet  problem  in  Q  =  (0,1)^  with  data  consistent  with  the  exact  solu¬ 
tion  Uj(xi,X3)  =  U2(xi,Z2)  =  sin()rxi)sin(TX3)  and  computed  its  finite  elements 
solution  on  a  9  X  9  uniform  grid  of  serendipty  squares.  For  this  finite  element 

solution  we  computed  the  values  of  the  relative  error  6(Xr>  a) 

(4  X  4)  Gauss-Legendre  points  and  at  the  f7%-superconvergence  points  given  in 
Table  3a  for  the  central  element  of  the  (9  x  9)  square  mesh  which  coincides  with 
/4  5\’ 

the  square  - j  .  In  Table  6a  we  give  the  values  of  the  relative  error  at  the 

(4  X  4)  Gauss-L^endre  points  while  in  Table  6b  we  give  the  values  of  the  error  at 
the  17%-superconvergence  points  given  in  Table  3a.  It  can  be  seen  that  the  rela¬ 
tive  error  at  several  of  the  (4  x  4)  Gauss-Legendre  points  is  nearly  47%  while  at 
the  f7%-superconvergence  points  from  Table  3a  the  relative  error  does  not  exceed 
3%.  Thus,  the  f7%-superconvergence  points  of  Table  3a  should  be  used  as  sam¬ 
pling  points  in  the  quartic  serendipity  element  instead  of  the  points  of  the  3x3 
Gauss-Legendre  product-rule. 


8  Summary  of  conclusions 

1.  We  presented  a  study  of  superconvergence  for  finite  element  approximations 
of  plane  elasticity.  We  employed  a  computer-based  methodology  which  takes 
directly  into  account  the  topology  of  the  grid,  the  element-space,  the  class  of 
solutions  and  the  value  of  Poisson’s  ratio. 

2.  We  determined  the  superconvergence  points  for  the  components  of  the  gradi¬ 
ent  of  the  displacement,  the  strain  and  the  stress.  We  observed  the  followring: 

a.  For  meshes  of  triangles  of  degree  p,  1  <  p  <  4. 
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(i)  For  some  mesh-p&ttems  (Regular,  Criss-Cross,  Chevron)  there  exist 
superconvergence  points  for  the  components  of  the  gradient  of  the 
displacement.  These  points  are  the  same  as  the  superconvergence 
points  for  the  components  of  the  gradient  for  finite  element  solutions 
of  Poisson’s  equation  given  in  [1].  The  location  of  the  points  does 
not  depend  on  the  value  of  Poisson’s  ratio. 

(ii)  There  are  no  superconvergence  points  for  any  of  the  stress-components 
(except  for  the  normal  stress-components  when  the  Poisson’s  ratio 
is  equal  to  zero)  or  the  shear  strain  in  any  of  the  mesh-pattem. 

(iii)  For  p  =  4  there  are  no  superconvergence  points  for  any  of  the  quan¬ 
tities  in  any  of  the  patterns: 

(iv)  Suitable  sampling  points  for  the  stresses  can  be  obtained  by  locating 
minimal  ij%-superconvergence  points  in  each  pattern  for  each  stress- 
components,  for  each  element-degree  p. 

b.  For  meshes  of  squares  of  degree  p,  1  <  p  <  4. 

(i)  For  elements  of  the  tensor-product  or  the  intermediate  family  for 
1  S  P  ^  4  and  elements  of  the  serendipity  family  for  p  =  1  and  2 
the  points  of  the  p  x  p  Gauss-Legendre  product-rule  are  supercon¬ 
vergence  points  simultaneously  for  the  components  of  the  gradient 
of  the  displacement,  the  strain  and  the  stress. 

(ii)  For  the  cubic  serendipity  square  there  exist  four  superconvergence 
points  and  one  superconvergence  line  the  components  of  the  gradi¬ 
ent  of  the  displacement.  For  the  normal  components  of  stress  (for 
non-zero  values  of  Poisson’s  ratio)  and  the  shear-strain  there  is  only 
one  superconvergence  point  located  at  the  center  of  the  element. 

(iii)  For  the  quartic  serendipity  square  there  are  no  0%-superconvergence 
points  for  any  of  the  solution  quantities.  However  it  is  possible  to 
locate  »;%-superconvergence  points,  for  the  components  of  the  gra¬ 
dient  and  the  stress,  for  small  values  of  if%  (i/%  <  2.5%).  The 
values  of  the  solution  quantities  are  much  more  accurate,  asymp¬ 
totically,  at  the  corresponding  f7%-superconvergence  points  than  the 
values  of  the  quantities  at  the  4x4  points  of  the  Gauss-Legendre 
product-rule. 
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Appendix 

“Harmonic”  monomials  of  degree  p  + 1  (1  <  p  <  4) 

Below  we  give  the  “harmonic’'  basis  monomials  of  degree  (p  + 1),  for  the  space 
defined  in  Section  6,  for  the  equations  of  plane  elasticity  and  1  <  p  <  4. 

a.  Quadratic  “harmonic"  monomials 

Q\^  01,03,63, 63)  =  0|  xj  +  <*2*2  + 


(xj, X2, 0|, O3, 63, —  6j  X|  +  63 X3  +  ^  ^1^2 

where 

A  +  2/X  A  +  /i^  ,  X  +  2n  ^  A  +  /i 

b.  Cubic  “harmonic”  monomials 

Ql^  (X|,  X2;  O3, 04,  64)  =  Oj  X|  +  O3  X3  '  J3  X|X3  +  O4  XjX3 


(®li  *21  ®3»  ^4)  —  6|  xf  +  63  X3  +  ^  ®J*2  ^4  *1*2 


where 


o,  =  - 


03  =  - 


^  o 

k  a  '* 

3(A  +  2p)  ^  3(A  +  2/i) 

^  3(A  +  2p)  ^  3(A  +  2p) 

A  +  2p  A  +  p 

3,.  3^ 

,  A  +  2p  ^  A  +  p 

3p  ^  3,. 

c.  Quariic  “harmonic”  monomials 

^2^"(*i>  *2;  «3»  06»  *5)  =  <*1*1  +  <*2  *2  +  ®3  *1*2  +  <*4  *1*2  +  <*5  *1*2 

^2^"(*l>  *2'.  «3»  *s)  =  *1*1  +  *2  *2  +  *3  *1*2  +  *4  *1*2  +  *5  *?*2 

where 
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*4 


3/1  3(A  +  2/i) 

2(A +  /*)“»  2(A  +  /i) 


«  -  /* 

4(A  +  2/i)  “  6(A  +  2/i) 

(f.  Quintic  harmonic”  monomiala 


Zfi  ,  3(A  +  2/i)^ 

2(A  +  /i)**'  2(A  +  /i)‘^ 

^  “  4(A  +  2/<)  ‘  6(A  +  2/i)  ^ 


(2|, Xjj  a^, a5, 64, 65)  —  a|  if  +  flj if  + 13 *1*3  +  ^4  *1*3  +  <*5 *1*2  *1*2 


(?3^  (ij,  I^;  04,  Ag,  64, 65)  —  6|  if  +  63 13  +  ^  ®1^2  ^4  ®1®2  ^  ®1®2  ^6  *1®2 


where 


a,  =  - 


a-  =  - 


=  - 


a,  =  - 


A-f /i 


6«- 


2(A  +  2/i)  *  2(A  +  2/i) 


A  +  2/1  A  +  /i 

~  ®4  > 


<*5 »  ^  — 


A  +  /I 


-  k 

®5”  o/\  .  0..\®S 


2/i  2/i 


5(A  +  2/i)  «  10(A  +  2/i) 


(A  +  2/i  (A  +  /i)  ^ 

Ok  —  — 0^  , 


6«  = 
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Minimal  values  of  i/%  for  the  stress-components 

Degree  of 

Chevron  pattern 

Criss-Cross  pattern 

the  elements 

WESM 

<Tjj 

<^11 

<r,2 

p  =  2 

18.17% 

5.75% 

11.79% 

12.19% 

5.53%(t,) 

30.50%(tj) 

28.75%(t,) 

24.51%(ra) 

p  =  3 

8.71% 

13.33% 

11.5% 

3.88% 

10.77%(t,) 

1.07%(rj) 

6.91%(t,) 

6.07%(r2) 

Table  1.  ri%-superconvergence  of  stress-components  in  the  meshes  of  triangles: 
Minimal  values  of  ri%  for  the  components  tn  the  Regular,  Chevron  and 

Criss-Cross  mesh-pattems. 


29 


Table  2.  Suptrconvtrmnce  joints  for  cubic  serendipity  square  elements:  Super- 
convergence  points  for  — .  Note  that  there  are  four  superconvergence  points 

vX|  vX| 

and  one  superconvergence  line  and  are  valid  also  for  the  class  of  general  solutions 
and  all  values  of  Poisson’s  ratio. 
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Quartic  serendipity  square  elements 


t7%-superconvergence  points  for 

OX\ 

Poiiit 

*a 

•7% 

-0.6364 

-0.5273 

5 

0.5273 

0.150 

6 

0.6364 

0.952 

7 

2.154 

8 

2.154 

Tkble  Sa.  ij%’Superconvergence  points  for  quartic  serendipity  square 
Sampling  points  for  with  q%  <  S.5%  for  Poisson’s  ratio  v  —  0.30. 


dements: 
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Quartic  serendipity  square  elements 

7%-superconver^nce  points  for  Cu 

Poiiit 

B 

*2 

•7% 

1 

-0.6727 

0.724 

2 

-0.4909 

0.499 

3 

0.4909 

0.499 

4 

0.6727 

0.724 

Tftble  Sb.  ti%-8upeTXMnveryence  points  for  quartic  serendipity  square  elements: 
Sampling  points  for  with  q%  <  0.75%  for  Poisson*s  ratio  v  =  O.SO. 
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Quartic  sercndi|Mty  Kiuare  elemenU 

17%-superconvergence  points  for  0^ 

Poiot 

mm 

*» 

17% 

1 

-0.5818 

1.334 

2 

0.815 

3 

0.5636 

0.0000 

0.815 

4 

0.5818 

1.334 

Ikble  Sc.  ri%-»uperoonvergenee  points  for  fuartie  strtniipity  square  eiements: 
Sampling  points  for  with  q%  <  1.5%  for  Pmsson*s  ratio  v  =  O.SO. 
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Rate  of  convergence  at  the  superconvergence  points 

Regular  pattern;  quadratic  triangles 

h 

5-3  max 

*“'€OolV^*i  ax,  y 

.25 

0.032184 

2.057 

.125 

0.003927 

2.011 

.0625 

0.000488 

1.999 

.03125 

0.000061 

1.996 

Ikble  4a.  Rate  of  convergence  at  the  superconvergence  points:  Values  of 

max  1(3^  ~  ~  ^^)(*^)|*  I^irichlet  prob- 

lem  in  n  =  (0,1)^  with  data  consistent  with  the  exact  solution,  u,(x,,X2)  = 
^(xiyX))  s  sin(xx,)8in(xx]),  Oq  :=  (0.25,0.75)^.  Grids  of  quadratic  triangular 

elements  in  the  Regular  pattern.  Note  that  the  values  of  at  the  superconvn* 

vX| 

genoe  points  are  superconvergent  with  rate  of  convergence  equal  to  3. 
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Quartic  Mrendipity  square  elements 


17%-superconvergeQce  for 


max 


.12053504E+00 

30.8570 

.70559307E-02 

28.9011 

.43019091E-03 

28.1930 

Table  5.  Rate  convergence  at  the  superconvergence  points:  Values  of 

max  *****  1(3^  ~  ^^)(*^)|*  Dirichlet  prob- 

I 'dxi  a*!-''  'I  *f'€0,l'^axi  dxi-''  ^  'I 

lem  in  n  =  (0,1)^  with  data  consistent  with  the  exact  solution,  U}(X|,X2)  = 
ttjCxijXj)  s  sin(xxi)sin(xx3),  Qq  (0.25,0.75)’.  Meshes  of  cubic  serendipity 

du-t 

squares.  Note  that  the  values  of  <7^  at  the  superconvergence  points  are  supercon- 
vergent  with  rate  equal  to  4. 


Relative  error  at  the  4x4  Gauss-Legendre  points 

P<rat 

*1 

*2 

-.861136311600 

-.861136311600 

30.5947 

-.861136311600 

.861136311600 

30.5948 

-.861136311600 

-.339981043600 

30.4170 

-.861136311600 

.339981043600 

30.4171 

5 

.861136311600 

-.861136311600 

30.5948 

6 

.861136311600 

.861136311600 

30.5947 

7 

.861136311600 

-.339981043600 

30.4171 

8 

.861136311600 

.339981043600 

30.4170 

9 

-.339981043600 

-.861136311600 

6.1581 

10 

-.339981043600 

.861136311600 

6.1580 

11 

-.339981043600 

-.339981043600 

46.5763 

12 

-.339981043600 

.339981043600 

46.5762 

13 

.339981043600 

-.861136311600 

6.1580 

14 

.339981043600 

.861136311600 

6.1581 

15 

.339981043600 

-.339981043600 

46.5762 

16 

.339981043600 

.339981043600 

46.5763 

Tbble  6a.  Values  of  the  relative  error  6(«. 


du-, 

•  ♦ 


u,U;^,^,t)  at  the  Ix-l  Gauss- 


Legendre  points:  Quartic  serendipity  square  element.  Dirichlet  problem  in  H  = 
(0,1)^  with  data  consistent  with  the  exact  solution  U|(X|,Z3)  =  U2(xi,X2)  = 
sin(xxj)sin(xx2).  The  domain  1)  was  discretized  using  a  9x9  uniform  mesL  of 
quartic  serendipity  elements.  The  relative  errors  are  reported  for  the  element  at 
the  center  at  the  center  of  the  mesh.  Note  that  the  values  of  the  relative  error  at 
some  of  the  Gauss-Legendre  points  exceeds  45%. 
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Relative  error  at  the  points  from  Table  Sa 

Point 

*2 

.5797 

.5797 

-.63640 

.0001 

-.5273 

2.8799 

5 

.5273 

2.8799 

6 

.63640 

.0001 

7 

.5797 

8 

.5797 

'Table  6b.  Values  of  the  relative  error  ©(as,;  Uf^^  A,  t)  at  the  r}%-supercon-‘ 

CrXj 

vergence  points  (from  Table  S):  Quartic  serendipity  square  element.  Dirichlet 
problem  in  (I  —  (0,1)^  with  data  consistent  with  the  exact  solution  U|(X|,X3)  = 
=  8in(«‘Xi)  sin(xx2).  The  domain  0  was  discretized  using  a  9x9  uniform 
mesh  of  quartic  serendipity  elements.  The  relative  errors  are  reported  for  the 
element  at  the  center  at  the  center  of  the  mesh.  Note  that  the  values  of  the 
relative  error  do  not  exceed  3%. 
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List  of  Figures 

Fig.  1.  An  example  of  a  locally  periodic  grid  of  squares:  An  adaptive  grid  with  a 
periodic  mesh-subdomain.  The  mesh  outside  the  subdomain  was  refined  in  order 
to  control  the  pollution>enor  in  the  interior  of  the  subdomain. 

Fig.  2  Periodic  meshes  of  triangles,  (a)  IRegular  pattern;  (b)  Chevron  pattern;  (c) 
Union-Jack  pattern;  (d)  Criss-Cross  pattern. 


Fig.  S.  Superconvergence  points  for  for  the  class  of  harmonic”  solutions  of 

the  equations  of  plane  elasticity:  Triangular  elements  in  the  Regular  pattern*..  The 
superconvergence  points  are  located  at  the  intersection  of  the  contours  {Qi^] 

f ,  T),  t  =  1, . . . ,  4.  In  the  Figures  the  0%-contours  of  the  error  corresponding  to  the 
various  monomials  were  drawn  with  different  thicknesses  and  the  element  bound¬ 
aries  were  drawn  with  dashed  lines.  The  0%-contours  and  the  superconvergence 
point?  are  given  for:  (a)  Linear  elements  (p  =  1);  (b)  Quadratic  elements  (p  =  2); 
(c)  Cubic  elements  (p  =  3).  Note  that  the  superconvergence  points  are  shown  in 
each  Figure  by  a  solid  circle  and  are  also  superconvergence  points  for  the  class  of 
general  solutions  and  for  all  values  of  Poisson’s  ratio. 


Fig.  4.  Superconvergence  points  for  for  the  class  of  ‘harmonic”  solutions 

of  the  equations  of  plane  elasticity:  Triangular  elements  in  the  Criss-Cross  pat¬ 
tern.  The  superconvergence  points  are  located  at  the  intersection  of  the  contours 
{Qi^;f,f),  i  =  1,...,4.  In  the  Figures  the  0%-contours  of  the  error  cor- 

responding  to  the  various  monomials  were  drawn  with  different  thicknesses  and 
the  element  boundaries  were  drawn  with  dashed  lines.  The  0%-contours  and  the 
superconvergence  points  are  given  for:  (a)  Linear  elements  (p  =  1);  (b)  Quadratic 
elements  (p  =  2);  (c)  Cubic  elements  (p  =  3).  Note  that  the  superconvergence 
points  are  shown  in  each  Fig.  by  a  solid  circle  and  are  also  superconvergence  points 
for  the  class  of  general  solutions  and  for  all  values  of  Poisson’s  ratio. 


Fig.  5.  q%-superconvergence  regions  for  for  the  class  of  “harmonic"  so- 

vX| 

lutions  of  the  equations  of  plane  elasticity:  Triangular  elements  in  the  Regular 
pattern.  The  regions  are  given  for:  (a)  Linear  elements  (p  =  1); 

(b)  Quadratic  elements  (p  =  2);  (c)  Cubic  elements  (p  =  3).  The  i/%-levels  5%, 
15%,  30%  (dark,  light,  lighter  gray)  were  employed. 


Fig.  6.  q%-superconvergence  regions  for  for  the  class  of  “harmonic”  solutions 

vXj 

of  the  equations  of  plane  elasticity:  Triangular  elerr^’nts  in  the  Chevron  pattern. 


The  regions  f ;  f)  are  given  for:  (a)  Linear  elements  (p  =  1)  (17%-levels: 

5%,  15%,  30%);  (b)  Quadratic  elements  (p=  2)  (7%-levels:  40%,  50%,  60%);  (c) 
Cubic  elements,  (p  =  3)  (17%-levels:  5%,  15%,  30%). 


dtl 

Fig.  7.  Ti%-superconvergence  regions  for  for  the  class  of  “harmonic”  solutions 

vX| 

of  the  equations  of  plane  elasticity:  Triangular  elements  in  the  Criss-Cross  pattern. 
The  regions  t;  f)  are  given  for:  (a)  Linear  elements  (p  =  1)  (f7%-levels: 

5%,  15%,  30%);  (b)  Quadratic  elements  (p  =  2)  (*7%-levels:  10%,  30%,  50%);  (c) 
Cubic  elements  (p  =  3)  (17%-levels:  5%,  15%,  30%). 


Fig.  8.  q  96-superconvergence  regions  for  and  class  of  “harmonic” 

solutions  of  the  equations  of  plane  elasticity:  Triangular  elements  in  the  Regular 
pattern,  (a)  ^J*(Q“^;r,!f),  p  =  2,  min  17%=  18.17%;  (b)  7i3|^(Q“^;f,f'), 
p  =  2,  min  17%  =  5.75%;  (c)  7lJ*((2“^;f,lf),  p  =  3,  min  17%  =  8.71%;  (d) 
7i;J(Q“^;f,r),  p  =  3,  min  17%  =*13.33%.  The  f7%-levels  10%,  30%,  60%  (dark, 
light,  lighter  gray)  were  employed. 


Fig.  0.  q%-supereonvergence  regions  for  (Tu  and  for  the  class  of  “harmonic” 
solutions  of  the  equations  of  plane  elasticity:  Triangular  elements  in  the  Chevron 
(waern.  (a)  p  =  2,  min  7%  =  11.79%;  (b) 

p  =  2.  min  7%  =  12.19%;  (c)  p  —  min  7%  =  11.50%;  (d) 

f),  p  =  3,  min  17%  =  3.88%.  The  T7%-level8  10%,  30%,  60%  (dark, 
light,  lighter  gray)  were  employed. 


Fig.  10.  q%-supera)nvergence  regions  for  <t,j  and  Ojj  for  the  class  of  “harmonic” 
solutions  of  the  equations  of  plane  elasticity:  Triangular  elements  in  the  Criss- 
Cross  pattern,  (a)  TIJ*  (C"^;  t,  T),  P  =  2,  inin  7%  =  5.53%,  min 7%  =  30.50%;  (b) 

=  2, inin7%  =  28.75%, min7%  =  24.51%;  (c) 7J2*(C'‘"";f,r), 

p  =  3,  min 7%  =  10.77%,  min7%  =  1,07%;  (d)  7Jj’^(C“^;f,T’),  p  =  3,  min 7%  = 

T,  Tj  12  T, 

6.91%,  min 7%  =  6.07%.  The  7%-level8  10%,  30%,  60%  (dark,  light,  lighter  gray) 

‘’’2 


were  employed. 


Fig.  11.  Superconvergence  points  for  the  class  of  “harmonic”  solutions  of  tiie 
equations  of  plane  elasticity:  Cubic  serendipity  square  elements.  The  supercon¬ 
vergence  points  are  located  at  the  intersection  of  the  contours  (Q,“^;^7T), 

t  =  1,...,4.  Superconvergence  points  for  (a)  -5-^;  (b)  -r-^;  (c)  <7,,;  (d)  c.j.  Note 

CrX|  0X2 

that  for  the  components  of  the  gradient  there  are  four  superconvergence  points 
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and  one  superconvergence  line.  For  the  components  of  stress  and  the  shear-strain 
there  is  only  one  superconvergence  point  located  at  the  center  of  the  element. 
The  superconvergence  points  are  shown  in  each  Figure  by  a  solid  circle  and  are 
also  superconvergence  points  for  the  class  of  general  solutions  and  for  all  values  of 
Poisson’s  ratio. 


Fig.  12.  ii%’8uperconvergence  rtgiona  for  the  class  of  “harmonic"  solutions  of 
the  equations  of  plane  elasticity:  Cubic  serendipity  square  elements.  The  regions 

•  (*^)  ^5  (^)  (®)  ^11 5  (**)  *^12  Poisson’s 

dvL  Qm  ^  ^ 

ratio  1/  =  0.30.  For  -3-^,  the  i;%-level8  5%,  10%,  25%  (dark,  light,  lighter 

dX|  0X2 

gray)  were  employed;  for  <Tii  the  f;%-levels  25%,  50%,  75%  were  employed. 


Fig.  IS.  rf%~superconvergence  regions  for  the  class  of  “harmonic"  solutions  of 
the  equations  of  plane  elasticity:  Quartic  serendipity  square  elements.  The  regions 

ffven  for  F(u) :  (a)  (b)  (c)  o,,;  (d)  £,2  for  Poisson’s 

ratio  1/  =  0.30.  The  t/%-levels  5%,  10%,  25%  (dark,  light,  lighter  gray)  were 
employed. 


Fig.  14.  Common  25%~superconvergence  regions  for  all  Poisson*s  ratios  v,  0  < 
1/  <  0.35,  for  the  class  of  “harmonic"  solutions  of  the  equations  of  plane  elasticity: 
Quartic  serendipity  square  elements.  The  regions  Q  T")  are 

0<i/<0.35 

given  for  F(ti) :  (a)  (b)  <7„;  (c)  £„. 
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Fig.  4b 
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Fig.  4c 
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Fig.  5a 
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Fig.  5b 
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Fig.  5c 
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Fig.  6c 
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Fig.  7c 
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Fig.  8a 
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Fig.  8b 
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Fig.  8c 
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Fig.  8d 
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Fig.  9a 
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Fig.  9b 
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Fig.  9c 
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The  Labonitoiy  for  Nnneriail  Analysis  is  an  integral  part  of  the  Institute  for  Physical 
Science  and  Technology  of  the  University  of  Mar^nd,  under  the  gener?«I  administration  of  the 
Director,  Institute  fw  Physical  Science  and  Tedinology.  It  has  the  foUov  g  goals: 

•  To  conduct  researdi  in  the  mathematical  theoiy  and  computational  implementation  of 
numerical  analysis  and  related  topics,  with  emphasis  on  the  numerk^  treatment  of 
linear  and  nonlinear  differential  equations  and  problems  in  linear  and  nonlinear  algebra. 

•  To  help  bridge  gaps  between  computational  directions  in  engineering,  physics,  etc.,  and 
those  in  the  mathematical  community. 

•  To  provide  a  limited  consulting  service  in  all  areas  of  numerical  mathematics  to  the 
University  as  a  whole,  and  also  to  government  agencies  and  industries  in  the  State  of 
Maryland  and  the  Washington  Metropolitan  area. 

•  To  assist  with  the  education  of  numerical  analysts,  especially  at  the  postdoctoral  level, 
in  conjunction  with  the  Interdisciplinary  Allied  Mathematics  Program  and  the 
programs  of  the  Mathematics  and  Computer  Sdence  Departments.  This  includes  active 
collaboration  with  government  agendes  such  as  the  National  Institute  of  Standards  and 
Technology. 

•  To  be  an  international  center  of  study  and  research  for  foreign  students  in  numerical 
mathematics  who  are  supported  foreign  governments  or  exchange  agendes 
(Fulbright,  etc.). 

Further  information  may  be  obtained  from  Professor  I.  Babudka, Chairman,  Laboratory  for 
Numerical  Analysis,  Institute  for  Physical  Sdence  and  Technology,  University  of  Maryland,  College 
Park,  Maryland  20742-2431. 


